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1 Introduction 



Recently there has been a great development in understanding the physics 
of black holes in the context of string theory Q . The entropy of black holes, 
which is proportional to the area of the event horizon, is explained in string 
theory by counting the number of microstates composed of D-branes and 
open strings with the ends attached on them 0. The Hawking radiation of 
nonextremal holes is also described in string picture as the process of creation 
of a closed string by the collision of two open strings on D-branes ||. The 
emission rate can be calculated in string picture |4j] and agrees precisely with 
the result of semiclassical calculation for the case of scalar emission, including 
the correct expression for the greybody factors || . The D-brane description 
of black holes is good in the near-extremal well as the extremal case, 

and indeed it is good beyond the range of the value of couplings in which 
the D-brane description is a valid approximation. This indicates that there 
should be a universal 'effective string theory' which describes the low-energy 
dynamics of black holes, whether or not the coupling is weak || 0. 

In this paper we undertake the semiclassical calculation of the emission 
rate for the case of massless spin-| particles^. The black holes we will con- 
sider are those in D = 5, iV = 8 supergravity, which are obtained by aligning 
D5-branes, D-strings and momenta in type IIB superstring theory compact- 
ified on T 5 . The resulting expression for emission rate agrees with the one 
predicted from the effective string theory on the worldsheet of the bound 
states of D-strings and D5-branes. We can see from this expression that the 
decay process is interpreted as the collision of a boson and a fermion on the 
worldsheet that creates a fermion in spacetime. 



2 For four-dimensional black holes, the emission rate has been calculated in [^) for 
massless neutral fermions that have no couplings with the gauge fields. We think that 
it is appropriate to consider fermions in supergravity theory and take into account the 
nonminimal coupling with the U(l) field strengths when we make a comparison with 
D-branc calculations. 
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2 The Dirac equation in D = 5, N = 8 super- 
gravity 

In this section we investigate the Dirac equation in D = 5, iV = 8 supergrav- 
ity. What we want to obtain is the precise form of the radial wave equation 



( 2.40|) and the conserved flux /. 

Let us briefly review the basic property of D = 5, iV = 8 supergravity. 
This theory contains 8 gravitini (ip A ), 27 U(l) vector fields (A™), 48 Dirac 
fields {x ABC ) and 42 scalars besides the graviton. This theory has global E e 
and local USp(8) symmetries. 8 gravitini and 48 Dirac fields belong to the 
irreducible representations 8 and 48 of USp(8), respectively. 27 vector fields 
belong to the representation 27 of Eq. 42 scalars parameterize the coset space 
E 6 /USp(8) and these degrees of freedom are contained in a single quantity, 
the 27-bein V/j B relating the representations 27 of E e and USp(8). Note 
that there is no dynamical gauge fields corresponding to the local USp(8). 
The U Sp(8) connection is made out of the 27-bein, like the spin connection 
which is made out of the vielbein. 

The spin-| fields in D = 5, iV = 8 supergravity obey the following Dirac 



equation fi~2|| : 



l D mX ABC + | \r mn Fi J n V}fx BC]D n DE - (Trace)^l = 0, (2.1) 
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-pm — pa/„-l\m fmn — -p[m-pnl JpU ofl aU 
1 = 1 ( e )a > 1 = 1 1 > t mn = 2c *[m^ n ] 

where D m is the SO (1,4) and [/Sp(8)-covariant derivative and Qab is the 
?75'p(8)-invariant tensor. The term (Tr&ce) ABC means the subtraction of the 
trace part so that the whole left-hand-side become traceless with respect to 
the indices (ABC). 



We shall solve the equation (\2.1\ ) in the nonrotating, nonextremal, triply- 
charged black hole background. The form of the Einstein metric is well- 
known i 



ds 2 = 
A = 


-X'ahdt 2 + A^ 

fl fb fn i 


[h- l dr 2 + r 2 dQ 2 {3) 


Si = 
h = 


1 H — | sinh 2 di 

i_ r A 


(i = l,5,n), 



} 



(2.2) 
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In the charged black hole background, there are excitations in the U(l) 
gauge fields and the scalars also. To know the form for these fields we need 
some knowledge about the relation between the fields in the E 6 symmetric 
theory and the toroidally-compactified type II theory. 

We first convert the USp(8) indices into the 5*0(5) indices. From the 
Type IIB on T 5 viewpoint, the USp(8) index originates from the spinor 
index of the compactified 5*0(5) times a factor two, the multiplicity of SUSY 
in ten dimensions. We denote the representation 8 of U Sp(8) as 

^ A = (r^a), (2-3) 

and at the same time we introduce the 50(5) T-matrixQ = (r,)^ and 
the charge conjugation matrix C = (C) a p, C~ l = (0 _1 ) a/3 satisfying the 
standard identities: 





= (r a ) f 


{r 4 ,rj 


= 25 dfe 


or.o- 1 


"pT p* 

a a 


cc ] 


= -CC* = 1 



(2.4) 

As can be easily verified, C and CTa are antisymmetric while CT ^ are sym- 
metric. The U Sp(8) invariant tensor Q AB is defined as 

_ / n a M _ / o -sP a \ , , 

ilAB ~ \ n<# ) ~ [ si o ) ■ (2 ' 5) 

We shall then decompose the irreducible representations of USp(8) into 
those of 50(5). 

Firstly, the representation 27 of USp(8) is the second-rank antisymmetric, 
pseudoreal traceless tensor F AB that satisfies 

{F AB y = n AC n BD F CD , 

n AB F AB = 0. (2.6) 



3 We use the dotted alphabets a, b, ■ ■ ■ and p,q, - ■ ■ for tangent and world indices of 
internal T 5 . The Greek indices are used for the spinor indices of internal 50(5). 
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The quantity F AB = V/j B F^ n transforms as 27 of U Sp(8) and we translate 
it to the quantity with 5*0(5) indices as 

pAB lpo.fi pa p j3 p \ 

F af 3mn = -| (jrnn{C)<±P + J m na( C T a ) a ^j , 

^mn = 4 {Fmn^aTp + ^ F ^n{^ ohT ■ ( 2 - 7 ) 



Similarly, the spin \ field x ABC belongs to 48 of USp(S). It is antisym- 
metric and traceless with respect to the indices (ABC) 

X ABC = X IAB C] ^ABC = Q? (2g) 

and obeys the symplectic Majorana condition 



x ABC = c ■ n AA ,n BB ,n C ax A 6 , (2.9) 

where the bar means the Dirac conjugation and C is the charge conjugation 
matrix acting on the external (SO (1,4) spinor index. Again we decompose 
them into the irreducible representations of SO (5) as follows: 

x = (x a ^' i X 7; X (3y> 

X aPl = {{C- l ) at5 {C- l y & - {YaC- 1 ) aP (Y il C- 1 )' 1& )r l5 

x a \ = (r.c- 1 )^ - (c- 1 )^^ 

Xa/3-y = {CafiC^S — {CTa)al3{CT a ) 1 ^ ^ S (2-10) 

where we used the fact that the quantity 

are completely antisymmetric in indices (a(3 , ~f5). The spinor fields £ aa ^ aa 
and C, a ,f] a are identified, from the IIB on T 5 viewpoint, as the two gravitini 
with the internal vector index and the two dilatini, respectively. 
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We shall then investigate the 27-bein V U B . The infinitesimal E$ trans- 
formation on the representation 27 which is not in its USp(8) subgroup is 
expressed as 

5F u = y} j kl F kl , (2.11) 
where S /J KL has the property 



IKL 

r 

^IJKL\* 



Z/ KL = 0, 



)* = Xijkl, (2.12) 

and the indices /, J,K,L,--- are raised and lowered by flu, Q IJ . Rewriting 
the components in terms of 5*0(5) tensors, they become 

£a/?7<5 = S |(C) Q , j a(C) 7 5 — (CTa) a i3(CT a )~ 1 s} 
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^ b \CT bc ) T& 



i—X\t5 



Kpj = {{CU{C) 1T -{CT,U{CT% T }-s h \T bt C- 

= it a {(c- 1 )^(cr«) 7j -(r«c- 1 ) ft ^) 7i } 

+h"\T {a C- 1 r^CT h) ) 1 s - hKC-'r^CU. (2.13) 

Note that the number of parameters is 42, the correct number of scalar fields 
in D = 5, iV = 8 supergravity. We can recognize these parameters as the 
scalar degrees of freedom in type IIB on T 5 , namely, s and s* as the dilaton 
((f)) and the axion (x), s ab and s* ab as the components of NS-NS and R-R 
2 forms with two internal indices^ (Bpq,B' p ^), t a as the components of self- 
dual 4 form with four internal indices (Cp qi ,s) and h ab as the fiinfbein of the 
internal five dimensions (E£). 

In the present situation the fields B pq , B'^ q and C^s are n °t excited and 
we concentrate on the s, s* and ^-transformations. After some calculation, 



4 We distinguish the R-R 2 form (and the quantities related to it) from the NS-NS ones 
by adding a prime'. 
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the explicit form of the 27-bein Vfj B is expressed in terms of the relation 

pAB _ yABrpIJ . 
ran V IJ mn' 

Ft n = E^E h .Fl n , (E = det(£|)) (2.14) 

F± = E$4f* (2.15) 

f J T T a ) = Eie W ( S *)( ^np + iHL^ \ (2 lg) 

\ J mna J V / V n ' mn V in mnp ) 



We identify the fields 



pt> pm H H' K K' 

L mm 1 mn) 11 mnp) 11 mnpi ly mn, ly mn 



with those in toroidally compactified Type IIB theory coupled to KK, D3, 
NS1, Dl, NS5 and D5 charges. Note that s parameterize SU(1,1)/U(1) 
rather than SL(2, R)/SO(2), and this is why {H mnp , H' mnj) ) and (K mn , K' mn ) 
appear in combinations in the above expressions. 

Rewriting the action of D = 5, iV = 8 supergravity 

C = R-lG IJ , KL F^ n F mnKL + ^-d m Gj, lKL d m G IJ ' KL + ..., (2.18) 



^ 2~""- mn 24 

Gu,kl = ^ab^cdVij B Vkl (2-19) 
in terms of the fields with SO (5) indices, and inserting 



i 

2 ' 



s = — H ■ — K =0 

° ? A± mnp - ti mn w ; 



the action becomes 

^L=£ = R(g) + \(d m G^d m G M + d m Gd m G- 1 )-^d m <j ) d m( j ) 
— - { CCt ■ ■ F p p mn, 3 _i_ P <t>r<m u' ff' mn 

+e-*GK' mn K' mn ) + • • • , (2.20) 
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On the contrary, ten-dimensional type IIB action can be written (if we 
drop the four form potential) asQ 

S = J d w xsf^^R{g) + -tr{d m M- x d m M) 

~(H H')M-'^, )+••]> (2-21) 

M = *[\ (2.22) 
This action can be brought to the same form as ( 2.20|) by applying the fol- 



lowing procedure: we first dimensionally reduce it to five dimensions, that 
is, ignore the x$„.$ dependence and insert 

(~ a Arrirpa \ 
e m m m \ (n nn\ 

o Ei J' (2 - 23) 

~9mn = e^e n b Vab , G pq = % & E 4 \ h , ip = ~ In det G pq . (2.24) 
We then dualize the three-form field strengths in five dimensions 



H \ 1 / K q 



r 



rji i . ) ^^— e ^ Me mnpqr I js- lqr I (2.25) 

mnp j 9 \ J 

K mn = 2d[ m C n ] (2.26) 

K' mn = 2d [m C' n] , (2.27) 



and Weyl-rescale the external and internal metrics 

2 ~ 2 

<?mn — C^Qmn j Gpq = e ^Gpq- (2.28) 

Inserting \ = K mn = H mnp = into the action thus obtained, we arrive at 
the same expression as ( |2.20|) . 

5 Here the hatted vector indices run 0, 1, • • • , 9, while the dotted and undotted indices 
run 5, • • • , 9 and 0, 1, • • • , 4, respectively. 
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From the correspondence explained above, we know the form of all the 
bosonic fields in Eq symmetric theory. The fields other than metric take the 
following form: 

«,-*<*-*») = /r l /s 
l i 

^55 — ■ ■ ■ — <J- gg — Ji J 5 Jn 

_ £ I 
^99 = /l 6 fi fn 

detG M = h l flfn\ (2.29) 



i 

' 3 

2 
In 



2 



9 Tq sinh a n cosh a n 

~~ r 2 _i_ r 2 

~r~ n 

, Tq sinh a.\ cosh cti 

09 ,y«2 j 2 

sinh a 5 cosh a 5 
C o - ^| • (2-30) 

Putting all these into (|2.1|) , we obtain the Dirac equation for the fields 

r-.e-jr-jf- 

It is rather easy to see that the USp(8) connection in the covariant derivative 
vanishes, but calculating the interaction term is quite tedious. After some 
calculation, we have found that the interaction term is diagonal only for a 
subset of the spinor fields. Define 

^\t ) = e a ±ir 1 ho '-- £ (r b r a )y^±^), (2.3i) 

4 a=5,-,8 

then the Dirac equation for these fields becomes 

±e*V^JCd(T*)>&- y/G^G^ n (T%^] = 0, (2.32) 
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or, denoting if)?±s generally as ip an d letting e% = ±1 for and e 2 the 
eigenvalue of r 9 , this equation can be rewritten as 



t m D m + -t mn {-e x e-iVGK, 



mn 



e 2 \ G 1 GqqF. 



9 

99 1 mn 



0. (2.33) 



Note that these fields are interpreted, from the type IIB on T 5 viewpoint, as 
the gravitini whose vector indices are those of internal T 4 . 

In the present case where three gauge fields are excited, there are three 
contributions in the interaction term. For the other subset of the spinor fields 
the three contributions cannot be diagonalized simultaneously and the Dirac 
equation for them will be much more difficult to solve. 

In the following we shall consider ( f2.33j ) with e\ = e 2 = — 1. Putting 
(|3), (gggD and ( g30p into flOg) , it becomes as 



- r\h-^r d ip = {rT% + ((h* - l)x% + C + iAT°) } Uip 
where the index i runs 1, • • • , 4 and C is a function of r and 



(2.34) 



A = ^ (/ 1 - 1 sinh2a 1 + / 5 - 1 sinh2« 5 + /- 1 sinh2« n 



n 



r 



(2.35) 
(2.36) 



Note that we can change the function C arbitrarily by rescaling ip. 

Let us see the rotational invariance of ( p. 34 ). Defining the orbital angular 
momentum operator and the spin operator as 



L 



3*j ®i ) 



-^[r.r,], 



(2.37) 



the Hamiltonian in (|2.34j) commutes with the total angular momentum J, 



ij 



Lij + Sy . Therefore we shall obtain the solution as an eigenfunction of J = 
Jij Jij and, say, J12, Jsa- The solution can certainly be expressed as a linear 
combination of eigenf unctions of {L 2 , J 2 , Ju, J34}, or, two spins (J,j') of 
SU(2) x SU(2) ~ 50(4) and their ^-components (m,m r ) besides the orbital 
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angular momentum I. We denote them as ^y mm /- The basic properties of 
^jfmm' are listed in the appendix. 



Rewriting the Dirac equation ( 2.34 ) with the operators (E • L) and IT, it 
becomes, after a suitable rescaling, as 



dr v J 2 



Putting the form of the general solution 



ij, = e~ tujt ■ !f 1±1 Ml , , + G l±1 Mti) , } 
into ( |2.38| ), we obtain the radial wave equation 



(2.38) 



(2.39) 



—ujr\^h~^Gi + i t i 



r 1 d ,, 3 . s 
rh* — + (/ + - + A) 
dr 2 

, i d 3 . . 

rA ,__ (i + _ +A) 



Gi+i,i 

^2+1,2- 



(2.40) 



We can derive the equations for F^+i, G^+i from the above equations by 
substituting {F l>l+1 , G^+i} for {—G l+l ih F l+ ljl } and -A for A. 

Finally, for spinor fields satisfying ( |2.38| ), the following equation holds: 

d (| dQ^xh-^To^j = -j- (| dfi (3 ) • 

So we define the conserved flux as / = / df2(3) , 0IT0. / is expressed by the 
radial wave functions i^+i.j, Gz+i^ as 



/ = i(G* l+l l F l+ i^ - F* +l l Gi +li 



(2.41) 



3 The solution 

Let us solve the radial wave equation ( |2.40| ) . As in the case of scalar field, the 
equation cannot be solved analytically and what we shall do in the following is 
to divide the range of r in two (or more) regions and obtain the approximate 
solution in each region and then continue them. We also limit ourselves to 
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the 'dilute gas approximation' or the case in which two of the three boost 
parameters, namely, ot\ and a^, are large. 

To begin with, we define some radial variables for convenience: 



^1 + cosh 2p) , y = tanh p , x = 1 — y , 



(3.1) 



and we recall the relation between the boost parameters and the mass and 
temperature of the black holefJJ: 



1 



1 



2nT R 
1 

2tfF L 
1 

2nT H 
1 

2ttT_ 
M 



2ttT h 2ttT_ 
1 1 



2nT H 2vrT_ 

r cosh a\ cosh «5 cosh a r 

r 1 sinh a.\ sinh sinh a r , 



(cosh 2«i + cosh 2a 5 + cosh 2a n ) 



(3.2) 



where the temperatures T H , T_ are defined by (^) times the surface acceler- 
ations of the outer and inner event horizons. In the dilute gas approximation 
the temperatures can be written as 



Tr i = !^ e KI+Kl. sinh | a I 



tttq 
2 

i _ ^ojail+Kl-Kl 



(3.3) 



In terms of the new variable the functions /i, A, A are expressed as 



h 
A 

A 



2nr T H 
1-y 2 d 



I — y 2 tanh 2 «i) (1 — y 2 tanh 2 «s) (1 — y 2 tanh 2 a 
l + y tanh «i 1 + y tanh a 5 1 + y tanh a n 



4 cfy 



In 



1 — y tanh ail — y tanh 0:5 1 — y tanh a r 
We define the new set of radial wave functions Xi, X 2 as 
X x = (V 2 /V^F l+hl , 



(3.4) 
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x 2 

Vi 

Vo 



{V 1 /V 2 )*G l+1 , 



(3.5) 



•(1 + ytanhai)(l + ytanha5)(l + ytanha n ) 
ZttIh 



■(1 — y tanh«i)(l — y tanha5)(l — ytanha n ). 
where q is an arbitrary constant. Then the equations for Xi,X 2 are 

X 2 = V X X X 



d I + q 



' l-y* 
d I + q 
dy l-y 2 



-V 2 X 2 . 



(3.6) 



3.1 The horizon region 

We first solve the equation in 'the horizon region' or the region of y in which 

| cothai i5 | — 1 <C (1 — y) 



or 



r <r sinh ai j5 . 
In this region it is a good approximation to set 

tanhax = sgn«i , tanha 5 = sgna 5 . 



Then the equation (3^) with 

q = 0(ai) + 0(a 5 ) + 0(a n ) 
becomes, under the dilute gas approximation, 



+ (/ + ,) 



X 2 



2/ -1 (l-^) y 

-y-\l + y^ 
sgna n 









2nT_J 


( " 


ysu 




2ixT_ 



W(3.7) 
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or, defining X± = X\ ± 2X2, the equation for them becomes 



d 


r - 




dy 


[27rT H y 


2ttT_ 


d + ^ 
dy 


f " 






2ttT_ 



*4 



l + q- 
l + q + 



2lfF L 
2nTr 



X+. (3.8) 



The equation ( |3.8|) can be solved without further approximation. The solu- 
tion with the appropriate boundary condition is as follows: 



X + = y 1 ^(l-y 2 )^F(a + l,A7+l;2/ 2 )V2 7 

x _ = y-^(i-yZ) l -¥ F (a,p,r,y 2 )V2, 



4vrTr 



l + q 



lit) 



l + q + l 



ILO 



1 



7 2 2itT h ' 



(3.9) 
(3.10) 



2 AtiT l ' ' 2 47rT fl 

The normalization of the solution is such that the absorbing flux at the 
horizon is / abs = 1. 

From the solution obtained above we can derive the leading order behavior 
of Fi +hh Gi + ij at small (1 - y). For a n > 0, 



1+1,1 



.Tj 



— 14 



'1 _-y2X-i(Z+?0+|). 

y / -pf 1+30+1 



p/1 _ tCJ - 
V 2 2ttT h - 



r(/ + g + i) 



-iui T*fl 



+(1-3/ 



r - 



2\|(Z4-go+|) 4ttT l j -V 2 27rT H 



47rTf, ' V 2 47rT fl - 



)r(-/-g -l) 



p/ — I— qo+t _ iui 
1 V 2 4vrT L 



p/ go _ iw ' 
\ 2 47rT fl - 



(3.11) 



G 



(1 _ y 2 )-K'+90-|) 



1 



4-ttTr ^2 2ttTh 



ir(/ + g ) 



p^ j+go+l 



itj Np/ l+go+2 _ jui 
4nT L ) L V 2 47rT fl 



?0 



+(i-y 

0(ai) + 0(a 5 ) 



2^G+<20+±). 



gx| - ^) r H - go) 

-j-go+l _ 



4vrT L 



ir(^o 



AttTr ■ 



, (3.12) 
(3.13) 



and the expression for a n < is obtained by Tl <-> Tr. 
String symmetry 
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Let us comment on the 'effective string symmetry' of ( |3.8| ). As was dis- 
cussed in |7j for the case of minimally coupled scalars, the approximate wave 
equation near the horizon should have a SL(2, R) x SL(2, R) symmetry if 
the dynamics and the microstructure of black holes are well described by 
the 'effective string' theory This also holds in the present case. Define a 
two-component function $(p, r, a) as 



$ = (0i,0 2 r, (3.i4) 

■z{r-(r)—r- 



h = e 5V ' 4 2 * t h 2 ^ T -(sinh2p)-2X_(p)^ 

j(T-<7)+f + 



>2 



e 2v "^4^^ + 2,1(^2^41^^ (3.15) 



and the SL(2,R) generators as 

= (3.16) 

1 • / -.9 1 , \ ( , 9 i d\ 
R 1 = -sm r + ff — + -cos(r + a) cothp— + tanhp— , 

2 op 2 \ or do J 

1 / x 9 1 • / s f , d , d\ 
R 2 = --cos(r + cr)— + -sm(r + cr) cothp— + tanhp— , 

2 op 2 \ Or do J 

-^1,2,3 = -Ri,2,3 with <t replaced by — a, (3.18) 
then the equation for $ can be written as 

2{L ■ £)$ = 2(-LiEi - L 2 S 2 + L 3 S 3 )$ = 1 + g ( Z + g ) $. (3.19) 
3.2 The asymptotic region 

Next we shall solve the equation ( p.40|) in the 'asymptotic region', where 



x — (1 — y) -C | cothcti^l — 1 

or 

r 3> r sinh ai 5. 



15 



In this region we shall treat ( [2.40Q with q = —I. ( |2.40| ) becomes as 



-x 2 

dy 2 



VlXi , ^-X 1 = -V 2 X 2 
dy 



(3.20) 



V'i 



Vo 



■2nT, 



H 



■(1 + ytanhai)(l + ytanha 5 )(l + ytanha r 



{^r)y-\i + y)- l -\i - y) 1 

•(1 — ytanhai)(l — ytanhasXl — ytanha n ). 
and the second-order differential equation can be easily derived: 



d 2 ( V' ' 

d? +ViV2+ \A 

dfi + VlV2 + [m 



2V 2 , 
2Vi 



(V 2 ~ 2 X 1 ) 



(vr 2 x 2 ) 



0. 



(3.21) 



(3.22) 



(3.23) 



where the prime represents the derivative with respect to y. Leaving only 
the leading two terms in the x-expansion, these become as 



dx 2 
d 2 



+ 



rluj 2 
x 3 



Mu 2 + 1 - (I + If 

Ax 2 



| r 2 u 2 | Muj 2 + 1 - (I + 2) 2 



{vpxs, 

(V^X 2 ) 



0. 
0. 



dx 2 x 3 Ax 2 
and the solution is expressed in terms of the Bessel functions as 

f = r u/V2x, 
Xi = vh' 1 (AJ l+1 ^ 2£{l+2) (0 + SJ_ (J+1)+26( , +2) (0) , 

X 2 = V^ 2 ^ 1 (AJi +2 -2e(1+1)(0 - -BJ_( Z+2 )+2e(2+l)(0) > 

Fi+u = (VrV^C 1 (AJ l+1 _ 2e{l+2) (0 + BJ_ (I+1)+2e( , +2) (£) 

G 1+1,1 = (ViV^) 3 ^" 1 (AJi + 2-2e(l+l)(€) - & J-(l+2)+2e(l+l)(€) 

Mlu 2 



(3.24) 
(3.25) 



4(Z + !)(/ + 2)' 



(3.26) 

(3.27) 
(3.28) 
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From the above solution we can calculate the incoming flux at infinity 



fir 



and the behavior of Gi+1,1 at small £ 

A2- 1 - 1 



\A- 



) l B\\ 



(3.29) 



Gu 



\r u}) 



r(/ + 2) 
r(/ + 3) 



.7,3 B2 l+1 , i 

£'+2 -| £~'~2 



r(-z) 

m l+2 

T(-l - 1 



(3.30) 



3.3 Matching the solutions 

We would like to connect the two solutions in the 'matching region' 



cothai , 



The first thing to do is to see whether the solution in the asymptotic region 
( |3.27| ) is good in the matching region also. Actually the differential equa- 
tions for X\ and X 2 in the asymptotic region (|3]22|) and flpD have regular 
singularities at zeroes or poles of V 2 and V\. Among others they have zeroes 
at x = 1 — coth 0:1,5 and x = 1 + coth 0:1,5 respectively and some of them are 
slightly below x = 0. We have to take into account the singularities arising 
from these zeroes near x = when we are in the matching region, though we 
could neglect them in the asymptotic region. The effect of these singularities 
is different according to the signs of oi 5. 



3.3.1 01 > , a 5 > 

In this case V\ does not have zeroes near x = and the differential equation 
for X 2 or ( |3.23| ) has no singularities near x — 0. Therefore the solutions for 
X 2 and Gi + i t i in the asymptotic region are also good in the matching region 
and we obtain the coefficients A, B by matching the solutions for G/+1,/ in 
the horizon region and the asymptotic region. They become, for a n > 0, as 

r(i-^)r(i + 2)r(i + 3) , fu ^ 



r(¥-^v)r(¥-i^) K J W 



2 4ttT l > V 2 4nT R > 

I -I) 

\irT R/ 



zT(i - -**L-)r(-l - 2)V(-l - V — 1 



B - JZ^'LKZ^r ~:' ^-Hr^(^) -(3.31) 



\-nT R > \ 2 4ttT l ■ 
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and Tr <-» Tl for a n < 0. On the contrary, the differential equation for X\ or 
(|3.22 ) has additional two regular singularities at small negative x besides an 
irregular singularity at x = 0. So the solutions for T/ +li / in the asymptotic 
region is no longer good in the matching region and actually there is no 
simple way to approximate Fi + n in the matching region. 



3.3.2 ax < , a 5 < 

This case is similar to the previous subsubsection. This time the solution for 
Fi + ij in the asymptotic region is good also in the matching region. Matching 
Fi +1)i gives, for a n > 0, 



.4 



B 



1 V 2 27rTff 



r(/ + i)r(/ + 2) 



r(i±i 



4ttT l ) L V 2 



2 l +l 



r u) 



-iY(l «*L_ 

61 V2 2ttT h 



1)T( 



4vrT fl 



r o; 



i+2 



4ttT l . 



(3.32) 



and Tr Tl for a n < 0. Approximating in the matching region 

is difficult because of the additional two singularities near x = in the 
differential equation for X 2 . 



3.3.3 at > , a 5 < 

In this case both Vi and V 2 have a zero near x = 0. Neither the solution for 
Fi+ij in the asymptotic region nor the one for Gi + ij are good in the matching 
region. So in the matching region we have to consider, instead of (|3.24|) and 
(3.25), the following differential equations: 



d 2 | r 2 uo 2 | Mm 2 + 1 - (/ + l) 2 



dx 2 8x 3 



Ax 2 
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4(x + 5i) 2 4x(x + 5i) 
^ r>2 Mcu 2 + 1 - (/ + 2f 
dx 2 8x 3 4x 2 

2(/ + 3) 



+ 



A(x + 5 5 ) 2 Ax(x + 5 5 ) 



(V 2 " 2 X 1 ) = 0, 



{V x ~ 2 X 2 ) 



(3.33) 
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<5 15 = | cothai j5 | — 1. (3.34) 



These equations are the same as ( ft.24 ) and ( |3.25D except for the restoration of 



the singularities near x = 0. Note that there is just one additional singularity 
for each of the above two differential equations. In the matching region we 
can neglect the second term in the L.H.S. of these equations and the resulting 
equations have two regular singularities at finite x. These are of analytically 
solvable form. The two independent solutions for i^+i,/ are 



i+i,i 



(TOK^i-Ci) -3 

■(-Gr a F(l + a - c, a, a - b + 1; Cf 1 ) 

(W 2 )ki f (i-Ci)^- 

f , > ,r(o-6 + i)r(i-c) 

+ (-Ci)^ c F(l + a - c, 1 + b - c, 2 - c; Ci 
T(a-b+l)T (c-l] 
T(a)T(c- 




(3.35) 



Ci = -x/8 x , a = -1-2 + e(2l + 3), b = e , c = -I + 2e(l + 2) , (3.36) 

and that obtained from the above expression by a <-» 6. The solutions for 
Crz+i,* are obtained from those for F t+ ij by Ci ~^ Cs = ~ x /^5 an d Z — > — Z — 3. 
Connecting the solutions in the horizon region, the matching region and the 
asymptotic region in order gives the value of A: for a n > 0, 

r(§-^-)r(/ + 2)r(/ + 2) 

V(1+2l — iuj \V( l +3 _ iu> \ 
1 v 2 4ttT l > \ 2 iirT R > 



A= "J E ^ (ro*)-'-* ( ~ ) , (3.37) 




and T R *-> T L for a n < 0. Matching Fj+i^ and C^+i,; give the same result 
for A in the dilute gas approximation, but not for B. But we do not have to 
take the mismatch of B very seriously because \B\ <C \A\ always holds. 



3.3.4 at < , a 5 > 

This case is similar to the previous subsubsection and we can obtain the 
result from the one in the previous case by 1 <-> 5. 
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3.4 The greybody factor 

From the results obtained in the previous subsections, we can calculate the 
greybody factor. Neglecting B being much smaller than A, it becomes as 



a 



abs 



n(l + !)(/ + 2) / abs _ vr 2 r 3 (/ + !)(/ + 2) 



\A\ 



r u 



(3.38) 



where the factor 7r (* +1 W +2 ) converts the partial-wave absorption probability 
to the absorption cross section of the partial wave. Under the dilute gas 
approximation, it can be written in a simple form: for a n > 0, 



abs 



2 2 3 -{ai+as+otn) 

7T r t 



(/ + !)(/ + 2) (r LU 



,2/ 



ilo \-p/ l+2+qo 
4ttT l ) L V 2 



r(l 



>r(/ + i + 



90+1 
2 



>r(/ + 2+ f ) 



_ 2-2i-2+2go 
2 

(3.39) 



2nT H < 

q = 6(ai) + 8(a 5 ), 

and Tr <-> T L for a n < 0. In this expression [x] is the Gauss symbol of x and 
[n] = [n + 1/2] = n for integer n. Recalling that the area of the outer event 
horizon is expressed as 

Ah = 27r 2 rQ cosh oi\ cosh a$ cosh a n , 

the first factor in the right hand side of ( |3.39| ) can be rewritten as 

2 -2 7r 2 r 3 e -(ai+a 5 +an) _ l/ige 

I _|_ e 2a„ 

and approaches to Ah when ai j5 are negative and a n is negatively large. If 
either or both of 0:1,5 are positive, the greybody factor becomes small and the 
emission from this fermion mode is suppressed. This is because the fermion 
mode we have been considering has negative D1,D5 and KK charges. 

In the following we shall consider the case of 0:1,5^ < 0. In this case the 
greybody factor becomes 



(2+1,0 
abs 



A 



H 



1 + e 



2a, , 



+ + 2) /r aA 2 ' 
4{/!(/ + l)!} 2 V 2 J 



AttTr) \ 2 



4ttT l . 



r( 



2-eTh' 



(3.40) 
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The form of cr a b s is strikingly analogous to the one for scalar field computed 
in |13|)[l6|. We expect that the same expression can be obtained from the 
effective string theory calculation as the collision process of R- and L- movers 
emitting a fermion into spacetime. The relevant interaction vertex must be 
a three-point coupling of spacetime fermion and one holomorphic and one 
antiholomorphic operators, each with weight ^jp and respectively. 

In the case of s-wave (1 = 0), 0^ becomes 

ff 0.o) - _AlL e ^ + 1 ^ ( 3 4 n 

abs l + e^ (e ^ + i)( e ^_i)4T L - [6 ^ L) 

Sending a n to a large negative value, it becomes as 

«Sf = (3-42) 



The decay rate is expressed in general 

rM (w) = E|4^- (3-43) 

and, in the s-wave case and for a large negative a n we obtain the Hawking 
formula by summing over polarizations. 

rV)^ = t^^t- (3-44) 



4 Conclusion 

We would like to conclude this paper by reviewing the calculation and sum- 
marizing the results. We considered the emission of spin-| particles from 
charged, nonrotating black holes in five-dimensional N = 8 supergravity. 
The Dirac fields in this theory have a 'Pauli term' or a nonminimal coupling 
with £7(1) field strengths and through this coupling the emitted Dirac parti- 
cles can carry the charges out from the black hole. This interaction term can 
be diagonalized for a certain subset of spinor fields which are interpreted, 
from the type IIB on T 5 viewpoint, as gravitini whose vector indices are 
those of internal T 4 . We derived the radial wave equation and solved it using 
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the similar method as in the case of scalar field. The wave equation turned 
out to have a SL(2,R) x SL(2,R) symmetry at the horizon region, as in 
the case of scalar field and this property is considered to be crucial for the 
'effective string' interpretation. The problem of matching the solutions was 
much more complicated than the scalar case. We finally obtained the ab- 
sorption cross section or the greybody factor. The expression is very similar 
to the one for scalar field and we expect that it can be re-derived from the 
effective string theory calculation. In the near-extremal limit the greybody 
factor approaches to the horizon area times an unusual factor 1/2. This is 
canceled by another factor 2 that arises from summing over spin degrees of 
freedom of the emitted fermion and the Hawking formula are obtained. 

Let us note that, as we said in the previous section, the enhancement 
or suppression of the emission according to the charges of the black hole is 
due to the charges of the emitted fermion. The fermion mode we have been 
considering is the fermion with negative KK,D1,D5 charges. Actually we 
have discarded the seven alternatives, namely, the fermion modes with the 
other signs of KK,D1,D5 charges by choosing e\ = €2 = —1 in ( |2.33[ ) and 



taking Gi+ij) rather than (i^j+i, Unless aji,5, n are all negative 

the emission from the mode we have considered are suppressed, but there 
emerges another mode whose emission are enhanced. 

We are not sure if the greybody factor ( |3.39| ) can be fully explained, 
including the overall numerical factor, by the effective string theory, though 
it is strongly expected. To see the full coincidence with the effective string 
theory calculation, we need more detailed knowledge about its worldsheet 
action, in particular its fermionic degrees of freedom and their interaction 
with spacetime fields. This will be our future problem. 
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A The properties of ^y mm / 

In this appendix we list some basic properties of 



^%mm' nas the following eigenvalues: 



jj'mm' 



J 2 

L 2 
E 2 

E • L 



The operator II 
value of I by ±1. 



4j(j + l)+4j'(j> + l), 

21(1 + 2) 

3, 

I 



= < 



f (0 



for 

Z-2 for *Sl, „*?L , 



(A.l) 



^(r • x) commutes with J^- and II 2 = 1. II shifts the 



(i±i) 



(0 



j t±i. 



(A.2) 



T° commutes with J {j and the Hamiltonian. (r ) 2 = —1 and we set its 
phase as 



2 2 

J- * ; (±i 



— ^ i±i j , 

2 2 

— l ^ I l±l mr 



(A.3) 



The orthonormality of ^fji mm t is set as 



dfy 3 )* t(0 

We finally list some identities: 

(r-x)(r-9) 



(r-9)(r.x) 
n(E • l) + (E • L)n = -3n. 



r-^ + 4 + (E-L) 

ar 



(A.4) 



(A.5) 

(A.6) 
(A.7) 
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B The relation between the plane wave and 
the spherical wave 

In this appendix we shall review the plane wave and the partial wave solu- 
tions of the Klein-Gordon and the Dirac equations in a flat five-dimensional 
spacetime. 

We begin with the Klein-Gordon equation. The plane wave solution is 
expressed as 

<j> k (x) = e 

and its orthonormality are expressed as 



ik-x 



(B.l) 



/ dx 4 <P* k (x) ■ 4> k ,{x) = (2tt)V(A; - k'). 
The spherical-wave solution can be written as 

4'[w,l,m\,m2\ \ x 



(B.2) 



(B.3) 



where are the four-dimensional spherical harmonics and the indices 

represent their eigenvalues of the orbital angular momentum: 



L L -Y (t) 

LJ ij 1J i] 1 mim 2 



21(1 + 2)r« 



i(L 12 + L 34 )ll i i ) m2 = rmY®^ 
I(L 12 -L 34 )y« = m2Y®. 



(B.4) 



Its orthonormality is expressed as 

/ dxA ^\w,i, mu m 2 ]i. x ) ( t ) W,v^ m ' 2 }{.x) = (2t:)5(uj - J)5i V 5 mimli 5 m ^ m ' 2 (B.5) 

We consider the Dirac equation next. We are going to work in the repre- 
sentation where 



/ 1 



2S 



-1 



12 



\ 



"I J 



2S 



34 



-1 



-1 7 



(B.6) 
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i 

V • • 



, r° = -r 



-i ■ ) 



V 



(B.7) 



-i I 



where E^ = — |[Tj,rj]. The normalized massless plane wave solutions are 
obtained as the four column vectors in the (4 x 4) matrix 



V2k° 



[tr m d m y k - x . 



(B.8) 



Note that only two of the four column vectors are linearly independent. 
Denoting the two independent solutions with the given momentum k as ipk,i 
and ipk,2, the ortho normality is expressed as 



J d 4 x^l s (x)^,Ax) = {2v) A 5\k - k')5 s 



(B.9) 



When the momentum is along the Xi-axis the two linearly independent solu- 
tions are written as 



V2 



f 1 \ 
1 

V • J 



-ik°(t-xi) 



V2 



( ■ \ 
1 

V i / 



g-ik°(t-xl) ( R10 ) 



where the indices ± represent the eigenvalue ±1 of 2E 34 . The spherical wave 
solution is written as 



-,—iiot 



V2 



LOT 



where ^% mm i are the eigenspinors of the total angular momentum and they 



—Utlt 



■Ji+i{vr)W L i ±1 -— J i+2{vr)v L i ±1 , 



oor 



uor 



(B.ll) 



■ jj'mm 
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are expressed, in our representation of T matrix, as 



* i+ii, 



21+2 m—-^,m' 



/ / H-2m+l y(0 \ 




/ i-2m+l y(0 
\ V 2/+2 z m +l m ' / 







2«+4 m,m'+ 



Z-2m'+2 v^(H-l) 





2Z+4 m,m'-± 



/ 





Z-2m'+l y(0 

7 

l+2m'+l y(i) 





\ 



22+2 ■* m,m'+i 



22+2 ■* m,m' — i 



/ / l-2m+2 y('+l) \ 
V 2Z+4 2 m _I jm ' 



i til mm' 








i+2m+2 y(i+l) I 
22+4 2 m +i, m ' / 



(B.12) 



(B.13) 



Finally we would like to calculate the geometric factor appearing in (|3.39|) . 
It is calculated as 

l/i±l,I,mm'| 2 

m,m' 

where the coefficients f i+i l mm' are determined through 

2 '2' 

We define the polar coordinates r, 6 1 , ^ 2 as 

(Xl,X 2 ,X3,X 4 ) = 

(r cos # cos , r cos sin <f\,r sin # cos , r sin sin if 2) 
In this coordinate system the spherical harmonics have the form 

Recalling the decomposition formula 



(B.14) 



(B.15) 



i8 = ^ n (n+l)J- 3 

<rt V ^ 



4tt y/2 



uor 



J n+1 (tor) ■ X {n) (cos 9) (B.16) 
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X<»>(cos0) = Sin ^ +1) ^ } 
v^sinfl 

JdQ (3) \X^\ 2 = 1 (B.17) 



we can calculate the coefficient f i+i i by 

2 ? 2 



yj(l + 2m + + l)y^ jm ,X(') (cos^ cos ^i) (B.18) 
_ iv /(Z + 2m' + 2)(Z + 2)ll i+ m ! ) ;,X (!+1) (cos cos Vl ) 

Taking the reality of X® and the fact that 

J dn i3) Y^X^(cos9cos^) = J dn (3) Y^_ ml X^(cos9cos^) 
into consideration, the calculation becomes very simple and we obtain 

E \f^L mm ,\ 2 = - 3 {l + l){l + 2). (B.19) 
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